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Abstract 

There is a direct group-theoretical connection between neutrino mixing and horizontal symmetry 
that can be established without any dynamical input. Such a connection is reviewed and expanded 
in this article. For certain symmetry groups G including and 5*4, it is shown that a generic 
U{1) X Q Higgs potential of a valon yields exactly the alignments dictated by the group-theoretic 
approach, but energy can now be used to discriminate different alignments. This mechanism 
possibly explains why starting from an A4 group, the tribimaximal mixing matrix with an enhanced 
54 symmetry is more preferable than the one without it. 
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I. INTRODUCTION 



After Chadwick discovered the neutron in 1932, Heisenberg came up with the 'isotopic 
spin' to distinguish it from the proton. This new quantum number corresponds to a global 
symmetry SU (2) which explains the similarity of neutrons and protons in strong interactions. 
Now that we have three generations of quarks and leptons, should there not be a new 
'family quantum number' to distinguish them, and a corresponding 'horizontal (or family) 
symmetry' to explain the similarity of their interactions? 

Unlike neutrons and protons which have nearly the same mass, the masses of fermions of 
different generations are vastly different. Moreover, they mix. Thus if horizontal symmetry 
is present, it has to be badly broken, a fact which makes the symmetry difficult to recognize 
unless enough remnants survive from the breaking to tell the story. We suggest that such 
remnants are indeed present and they are hidden in the mixing matrix. 

Neutrino mixing is successfully described by the tribimaximal mixing (TBM) matrix [1]. 
This mixing is most frequently explained by models with an A4 symmetry, but models based 
on 5*4 as well as many other groups Q [2] can also do the job. With so many successful models, 
it is hard to know which of them is the correct horizontal symmetry, and indeed, whether 
horizontal symmetry really exists or not. 

To shed some light on this question from a different angle, we have previously developed 
a purely group-theoretical method to connect neutrino mixing with horizontal symmetry 
[3]. This method differs from the conventional one in that only symmetry considerations 
are involved. The presence of Higgs fields is never assumed, nor the details of a model 
Lagrangian. In this approach, the broken symmetry (which we call residual symmetry) and 
the minimal unbroken horizontal symmetry of the left-handed leptons are derived from the 
mixing matrix U. In particular, if neutrino mixing is given by the TBM, then the smallest 
unbroken horizontal symmetry is S4, not A4. Models based on can explain TBM because 
in those models there is an accidental symmetry which elevates the A4 symmetry to an 5*4 
symmetry. These results will be reviewed in Sec. II. 

The same formalism can in principle be applied to fermion mixing as well [3] , but unless 
we make suitable approximations to the CKM mixing, all that we get is the group SU{3), 
nothing smaller. For that reason we will not discuss quark mixing any further in this paper. 
This difference, however, brings up a very important question, as to whether there is a 
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common origin between neutrino mixing and quark mixing that is smaller than SU{3). This 
question will be taken up in a forthcoming paper [4]. 

The group-theoretical procedure can be reversed to obtain mixing matrices U from a 
symmetry group Q. The detail will be discussed in Sec. Ill, but let us outline now how 
that differs from the usual dynamical method. In the dynamical approach, U is to some 
extent determined by the vacuum alignments, and the alignments are in principle obtained 
by minimizing a Higgs potential V. There are however frequently too many ^-invariant 
terms available to construct V, so that just about any vacuum alignment can be obtained 
by a suitable choice of V. What one does then is to use the experimental results to decide 
what U one wants, and what alignments that can give rise to such an U, then proceeds to 
design a V that does the job. To keep it robust, this is often accomplished by imposing 
another symmetry in additional to Q to throw away the unwanted terms in V. 

In the group-theoretical approach, U is determined by the residual symmetries F (for the 
charged-lepton sector) and (i = 1, 2, 3) (for the neutrino sector) which are members of Q. 
No peeking of data is necessary so one is not prejudiced by the experimental outcome. There 
is no dynamics to tune, and the existence of Higgs is not assumed. For a given Q, there is 
usually more than one possible choice of residual symmetries and thus more than one U, 
but within equivalence only very few survives. By equivalence I mean the following. Since 
Major ana phases are not measurable in neutrino oscillation experiments, and since there is 
no way for group theory to distinguish which flavor neutrino is which, and which neutrino 
mass eigenstates are to be labeled 1, 2, or 3, we will consider two mixing matrices equivalent 
if they differ only by row and column phases, and by row and column permutations. With 
a convenient choice of permutation, the mixing can be specified by three mixing angles and 
one CP phase, contained in the triplet (sin ^12, sin ^23, sin ^ise""^). Two C/'s with the same 
triplet parameters are considered to be equivalent. 

Normally we need a F and three mutually commuting Gi to specify aU. A mixing matrix 
determined this way will be called a full-mixing matrix. These are the mixing matrices 
obtained from Q by symmetry alone. In order to allow TBM to occur under A4, we will also 
discuss partial-mixing matrices, which arc obtained from a full-mixing matrix by replacing 
two of its columns with arbitrary parameters, subject however to the constraint of unitarity. 
Equivalently, they are mixing matrices determined by F and one of the three Gj's, plus 
additional parameters needed to specify U. For each full-mixing there is always associated 
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three types of partial-mixings, which will normally not be mentioned separately. 

However, the mixing matrix obtained from a partial-mixing of Q with the other two 
columns filled has its own symmetry group Q which contains ^ as a subgroup. Unless the 
entries in the other two columns are carefully chosen, most likely Q would end up to be an 
infinite group, perhaps as large as SU{3). Thus, the entries of the other two columns are 
severely restricted if we stipulate ^ to be a finite group. In particular, if the entries are 
those given by the full- mixing matrix, then Q — Q. In the case of A4, the smallest finite 
group containing it is S^, and there is one partial- mixing mode of whose ^ is 5*4. That 
is how TBM is produced from A4,. 

In the case of 6*4, there arc two full-mixing patterns. One is the familiar TBM specified 
by the physical parameters (^, 0), and the other is one specified by 0), which is 

seldom mentioned in the literature. The full-mixing of is the Cabibbo- Wolfenstein mixing 
matrix, specified by (;^, ^-^5-), though TBM is allowed in a partial-mixing mode. There 
are two possible full-mixings in A5, but we will work out only one whose physical parameters 
are (:^, 0), where </? = (l-|--\/5)/2 = 1.618 is the golden ratio. Once again, these results 
are obtained in a purely group-theoretical way, without the help of any dynamics. 

Sees. IV, V, VI are devoted to the connection between group-theoretical and dynamical 
approaches. Defining an invariant eigenvector io be an eigenvector with eigenvalue -|-1, the 
basic connection is that the vacuum alignment in every IR has to be the invariant eigenvector 
of the corresponding residual symmetry operator in that IR [3] . That connection is reviewed 
in Sec. IV. In case such an eigenvector does not exist, the alignment would be forced to be 
zero and the coupling with such a valon cannot contribute. This in principle would reduce 
our ability to make models because we need at least three independent Yukawa couplings in 
each sector to fit the three charged-lepton and the three neutrino masses. Nevertheless, in 
the case of 5*4 and A^ which we discuss in some detail, there remain just enough parameters 
to fit all the data in type-II seesaw models. For A4 we also discuss the partial-mixing results, 
both because the phenomenologically interesting TBM belongs to that category, and because 
the result of the discussion would be useful in Sec. V. We see from these discussions that a 
large number of models can be constructed for each group, and indeed, many of them have 
already appeared in the literature. Similar discussions can be carried out for other groups. 

With the basic connection, one important question to ask is whether a generic Higgs 
potential can be constructed to yield just these vacuum alignments. As remarked before. 
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in order to avoid the appearance of arbitrary alignments, an additional symmetry would 
have to be imposed to throw away unwanted ^-invariant terms. In the case oi Q — A4, that 
additional symmetry turns out to be (an appropriate Z„ would do as well) as we shall 
show in Sec. V. The calculation involves writing down a generic U{1) x A4 potential in terms 
of CG coefficients in a convenient basis, then minimizing it. 

Can this result be generalized to other groups? On the one hand, this calculation requires 
intimate knowledge of the group A4, such as the exphcit CG coefficients and the details of 
all the X A4)-invariant terms, so it is difficult to see how to generahze it from A4 to 

other groups Q. On the other hand, the result of minimization can be expressed in general 
group-theoretic terms, as the invariant eigenvectors of members of A4, and that gives hope 
that one might find a way to bypass the details of to generalize the result to other groups 
Q. This is indeed true as will be discussed in Sec. VI. In the first part of that section, an 
alternative proof for the A4 result is presented that is quite independent of the specific A4 
details. This proof is then generahzed to other groups Q in the rest of that section. 

Besides demonstrating the consistency between the group-theoretical and the dynamical 
approaches, the introduction of the Higgs potential also serves another purpose. There is no 
way for group theory to choose between different residual symmetries, but Higgs potentials 
might, as different vacuum alignments might give rise to different Higgs energy. In that case 
the one with the lowest Higgs energy wins out. 

This mechanism can for example be used to solve a puzzle. TBM is a full-mixing of 5*4 
but only a partial- mixing of A4. The full-mixing of A4 is given by the Cabibbo-Wolfenstein 
matrix which gives too large a solar angle and too large a reactor angle. If A4 is the correct 
symmetry group, as many models seem to prefer, one has to answer two questions. (1) 
Why does Nature choose a partial mixing rather than a full-mixing more appropriate to 
A4 from the symmetry (group-theoretical) point of view? (2) Why should Nature choose a 
partial-mixing of A4 and then promotes it to a larger ^'4 to give TBM, rather than choosing 
5*4 directly as the symmetry group from the very start? 

With the result of Sec. V, we may have an answer to question (1). If the Higgs self 
couplings are all positive, then the full-mixing pattern has a higher energy than the partial- 
mixing pattern, which is presumably why Nature prefers partial mixing. Once in a partial- 
mixing mode, economy would like to make the mixing matrix TBM because ^ = 5*4 is the 
smallest group that contains A4. Question (2) is harder to answer, but the answer might 
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come from dynamics as well. In a forthcoming publication [4], it will be shown that if we 
ignore the much smaller quark-mixing with its third generation, then there is something 
common between quark mixing and leptonic mixing described by A4, but there is no such 
commonality between quark mixing and leptonic mixing described by 5*4. 

II. FROM MIXING TO HORIZONTAL GROUP 

We review in this section how the horizontal symmetry group Ql for left-handed leptons 
can be derived from the neutrino mixing matrix U [3]. We will assume to be a subgroup 
of SU{3), and we are particularly interested in small to minimize the amount of break- 
ing necessary to obtain the desired mixing. For that reason we will concentrate on finite 

subgroups of SU{3). 

Let Me be the 3x3 charged- lepton mass matrix, and Mg := M^Me the effective mass 
(squared) matrix connecting left-handed to left-handed charged leptons. This matrix is 
clearly hermitian. Similarly, let Mj, = be the Majorana mass matrix for the left-handed 
neutrinos. This matrix is symmetric on account of the Majorana nature of the neutrinos. 
In the basis where Mg is diagonal, the neutrino mixing matrix U is the matrix that renders 
U^MM diagonal. 

A residual- symmetry operator in the charged-lepton sector is a unitary matrix F which 
transforms Mg into itself: F^MgF — Me- Similarly, a residual-symmetry operator in the 
neutrino sector is a unitary matrix G which transforms M^, into itself: G^M,^G — M^,. G^ 
rather than G^ has to be used to maintain — mJ. An important consequence that 
follows is = 1. As a result, the eigenvalues of G are -1-1 or —1. Since detC — 1, G has 
one +1 and two —1 eigenvalues. There are exactly three such operators, given by 

Gi — Uiu[ — U2U2 — U^u\^ G2 — —Uiu[ -\- U2U2 — U^u\^ G3 — —UiUi — U2U2 + Usul, (1) 

where itj is the ith column vector of the mixing matrix U. It is easy to see that the three 
Gj's mutually commute, and that the product of any two is equal to the third. Moreover, if 
we define an invariant eigenvector to be an eigenvector with eigenvalue +1, then Ui is the 
invariant eigenvector of Gj. 

Since Me is diagonal, F may be taken to be any diagonal unitary matrix, with determinant 
1. If F is to be a member of a finite group, then there is an integer n so that = 1. We 
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divide F into two categories, the degenerate ones, and the non- degenerate ones. The latter 
consists of F's whose three diagonal entries are all different, which necessitates n > 3. This 
category is specially important because in that case Me must be diagonal whenever F is. 
This will be used in the next section to recover U from Ql. 

F and Gi are symmetry operators of the left-handed mass matrices, so they are the 
remaining symmetry operators of the left-handed leptons after the horizontal group Ql for 
left-handed leptons is broken. Conversely, the smallest Ql must be the group generated by 
F and G^, and this is how the group Ql can be obtained from the mixing matrix U. We 
shall use the notation = {F, G\, G2, G3} to denote this. Since the product of two distinct 
Gi is equal to the third, wc can write Ql = {F, Gi, Gj}, as long as i ^ j. We shall call the 
group generated this way the full group. 

The role of a non-degenerate F is to ensure that Mg is diagonal when F is, and the role 
of Gi is to use their invariant eigenvectors to construct the columns of the mixing matrix U. 
Since there is no need for more than one F to nail down Mg to its diagonal form, and since 
we want a minimal ^l, we include only one F in the generators of Ql- However, different 
choices of F will generally lead to different groups Ql- 

Sometimes we may be more sure of one column of U to be correct than another column. 
In that case we simply fill the other two columns of U with adjustable parameters, subject to 
the unitarity condition of U, of course. The full group Q associated with U then contains a 
subgroup Ql generated by {F, Gj}, where i denotes that particular column of U that we have 
confidence in. The group Ql will be called a partial group. There are clearly three possible 
partial groups, associated with i = 1,2, and 3 respectively. If for the sake of economy 
we demand ^ to be a finite group, then the entries of the other two columns are severely 
restricted. In particular, if Q = Ql, then the partial group is a full group. 

Note that since each Gi is of order 2, a partial group must have an even order, and a full 
group must have an order divisible by 4. Groups of odd order, such as A(27), would violate 
the assumption that at least one Gi is a member of the unbroken horizontal group. 

The symmetry of the right-handed leptons cannot be determined by this procedure, 
because mixing has nothing to do with right-handed fermions. Whenever it calls for a 
knowledge of those, we shall make the minimal assumption that they are invariant under Ql 
as well. In other words, if necessary, for simplicity we will implicitly assume Q — Ql this 
paper. 
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For tribimaximal mixing [1] , 



/ 2 a/2 \ 



V6 



-1 ^2 73 
V-1 V2 -y/3j 
the explicit form of Gi using (1) works out to be 



^0, 



(2) 



/ 1 



G 



10 



-2\ 



-2 -2 1 
V-2 1 -2/ 



/-I 



20 





n 





0\ 


; G30 — — 








1 




lo 


1 


0/ 



(3) 



2 -1 

V 2 2 - 

To obtain a minimal group, we shall choose for F the simplest non-degenerate matrix, 
Fq — diag(l,a;,a;^), where u — e}'^'^!^ . In that case, straight forward calculation using (3) 
shows that the full group is Qi, — {Fq, Gio, 6*20, G30} = 5'4, and the three partial groups are 

{-Po, Gio} = 5*4, {Fo, G20} = ^4, {Fo, G30} = 'S'a. 

Furthermore, it can be shown that as long as F is non-degenerate, no matter what n is 
(F" = 1), the full group Ql always contain 5*4 as a subgroup [3]. 



III. FROM HORIZONTAL GROUP TO MIXING 

Given a finite group Ql containing a non-degenerate element F, we can reverse the argu- 
ment of the last section to obtain all possible mixings associated with this group. Everything 
can be carried out in a purely group-theoretical manner, without the presence of Higgs fields, 
nor the help of a Lagrangian. 

In doing so, we may assume to be a full group, or a partial group. We will refer 
to the mixing matrices obtained with the former as full-mixing matrices of Ql, and those 
with the latter as partial-mixing matrices. Since partial-mixing matrices are obtained from 
full-mixing matrices by replacing two columns with arbitrary entries, subject only to the 
unitarity constraint, there is no need to discuss them separately. 

Another reason not to discuss them separately is that a partial-mixing matrix of Ql is a 
full-mixing matrix of a larger group Q, as mentioned in the Introduction. Hence we need to 
study only full-mixing matrices. 

The general procedure of constructing U from is as follows. First, separate the order-2 
elements in from those of order > 3. The former are candidates of G, and the latter are 
candidates of F. Next, identify other order-2 elements G' that commute with a G. Then 
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G" — GG' is an order-2 element that commutes with both G and G', so the normahzed in- 
variant eigenvectors of G, G' , G" in the F-diagonal basis constitute the three columns of the 
full-mixing matrix U. Which eigenvector occupies which column is a matter of convention, 
corresponding to different ways of labeling the three neutrino mass cigcnstatcs, or equiva- 
lently different ways of assigning G, G', G" to be Gi, G2, and G3. Moreover, the entries of the 
diagonal F may be permuted, resulting in permutations of the rows of U. In other words, 
the full-mixing matrix U can be determined only up to permutation of rows and columns, 
because group theory has no way of knowing which flavor neutrino is which, nor which mass 
neutrino eigenstate should be called 1, 2, or 3. When we present a U below, it is understood 
that such an ambiguity always exists. 

We must check that F, G, G', G" generate the group Qli and that the three normalized 
invariant eigenvectors in the F-diagonal basis are mutually orthogonal so that U is unitary. 

If we cannot flnd two order-2 elements in that mutually commute, which would be 
the case for example when the order of Ql is not divisible by 4, then Ql contains only 
partial-mixing matrices. The actual symmetry group Q of the resulting mixing matrix must 
be larger, with an order divisible by four. 

Since Majorana phases are not measurable in neutrino oscillation experiments, we will 
regard two C/'s differing only by Majorana phases to be equivalent. With that we may 
classify neutrino mixing also by the familiar Chau-Keung parametrization used in quark 
mixing, which by convention leaves the (11), (12), (23), (33) matrix elements of U real, 
through a suitable choice of the phases of the three columns and the three rows. A mixing 
matrix parametrized this way can be summarized by the triplet (sin 6'i2, sin ^23, sin 6*136"*''), 
exhibiting the mixing angles 9ij and the CP-phase 5. For TBM, this physical triplet is 

Although straight forward, the task of obtaining all the full-mixing matrices from a given 
group might seem rather daunting because of the large number of combinatorial choices of 
F and G's. In reality it is not as bad because mixings produced by members in the same 
conjugacy classes are often equivalent, i.e., they differ only by row and column phases, or 
the re-shuffling of rows and columns, so at least for relatively small groups Ql, there are 
only very few inequivalent mixings. 

With such a foresight, we should arrange group elements into conjugacy classes. One way 
to do so is to make use of the character table, although characters other than those in the 
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defining three-dimensional representation will not be used in this section. To illustrate the 
procedure, we shall discuss in the rest of this section how to obtain the full-mixing matrices 
from the groups — S4,A4^, and A^. These are the rotational symmetry groups of the 
regular polyhedrons. 

The elements of each of these three groups can all be obtained by the repeated multipli- 
cation of two members, a and b, which are known as the generators. In each case, a obeys 
a'^ — 1 and b obeys b^ — 1. These three groups differ from one another in the behavior of ab, 
which is of order 3, 4, 5 respectively for ^4, S4, and A5. In short, the presentation of these 
three groups are 

^4 = {a, b\a^ = 6^ = (ab)^ = 1), ,^4 = (a, b\a^ = 6^ = (ab)^ = 1), 

As = (a, b\a^ = b^ = {abf = 1). (4) 

A. S4 

5*4 is the group of permutation of four objects, with a character table given in Table 1. 
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Table 1. Character table of S4 



It contains five classes Ci (1 = 1, 2, • • • , 5) with \Ci\ elements, defined by the cycle structure 
of permutations shown in the second row. The order of the elements in each class can be 

read off from the cycle structure to be 1, 2, 2, 3, 4, respectively. 

There are five irreducible representations, x"; with dimensions 1, 1, 2, 3, 3. The two 
3-dimensional irreducible representations IR4 and IRS are very closely related. If g^' and 
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are the representations ol g & S^, respectively in IR4 and IRS, then g^ — g^' for g E C3 
or C4, and g^ — —g^' for g E C2 or C5. The defining three-dimensional representation (the 
matrices defining F and Gj) is IRS, not IR4, because = — 1 but = +1 for the class C2. 
Since G is of order 2, it belongs to €2- It also has eigenvalues +1, —1,-1, hence the trace of 
G, which is the character of C2, is —1 and not +1, so G must belong to IRS. 

The last row lists where the generators a and b belong. From the order of the classes, 
clearly b has to belong to C3, but a priori a could be in C2 or C3. C3 is excluded because then 
both a and b are even permutations, so the group they generate is the subgroup A4, not the 
whole S4. Hence a must be in C2. 

Since the naming of the four objects under permutation is arbitrary, without loss of 
generality we may take b = (134). To be a generator of 5*4 together with b, a must then 
contain the number '2', so there are three possibilities: ai=(12), 02 = (23) = baib^, and 
03 = (24) = b'^aib. For each Oj, one can find an a\ and an a'l so that members of the 
triplet Ai :— [oj, a^, a"] mutually commute and the product of two is the third. The explicit 
expression for the three triplets are 

Ai = [(12),(34),(12)(34)], = [(23), (14), (23)(14)], ^3 = [(24), (13), (24)(13)]. (S) 

Note that a^, e C2, but a" e C3. 

Next, let us see how to represent the generators as SU{3) matrices, in the basis where 
b — (134) is diagonal. Let us start off with b. Since b^ — 1 and its character Tr(6)=Tr(C4) 
vanishes, its diagonal form must be 6 = diag(l, cu, cu^) up to possible permutation of the 
entries, where cu = e^'^''-!'^ , This is the same as Fq in TBM. From (3) and the discussion 
below it, we know that Gio and Fq together generate 5*4, so we may take Gio to represent 
any one of the generators Oj. For definiteness we let it be a\ — (12). 

Please note that notations such as a, 6 may stand for the abstract elements defined by 
the permutation cycles, or matrices in various irreducible representations. 

To find the representation of the triplets Ai in the 6-diagonal basis, let B := (ai6)^ = 
(14)(23). Then BaiB^^ = which works out to be a'^ = G^q of (3). Furthermore, 
a'l = aia'i = G20 of (3). Hence Ai = [Gio, Gso, ^20]- Moreover, it is easy to check that 
A2 = Mift^ and A3 = b-'A.b. Hence A2 = Fo[G'io, G30, G2o]Fi and A3 = F^[Gio, ^30, G'2o]Fo- 

Any one of the triplets Ai is a possible choice ior G :— [G , G' , G"]. As to F, since F is of 
order > 3, it is either in C4 or C5. Let us consider these two cases separately. 
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1. F E . In cycle notation, we can write F — (xyz) for some numbers x, y, z chosen 
among 1,2,3,4. The permutation g = i ^ A brings F into diagonal form because 
gFg^^ = (134) = b = Fq. Since Fq together with any Ai generate 5*4, so does F and 
Af :— g~^Aig. The choice of G = {G, G', G") would then be any of the three . 

To compute U, we need to go to the basis where F is diagonal, in which case the 
representation of G is A^ in the 6-diagonal basis. We may therefore assume from the 
outset that F — b and G — Ai for some i. 

A similar step should be carried out for the cases to be considered later, but we will 
often skip that step and assume F directly to be an appropriate generator of the group. 

(a) G — Ai . If we let G — [Gi, G3, G2], then Gi — Gio and the full- mixing matrix is 
the TBM matrix Uq in (2). 

(b) G = A2 . If we let G = [Gi, Gs, Gs], then Gi = FoG^F^ and the full-mixing 
matrix is just the TBM matrix Uq, multiplied respectively by 1,0;, cu^ in the first, 
second, and third rows. Since row phases are adjustable, this is equivalent to the 
TBM in (2). 

(c) G = A3 . If we let G = [^1,^3,^2], then Gi = F^G^qFq and the full-mixing 
matrix is just the TBM, multiplied respectively by l,(jj^,tj in the first, second, 
and third rows. Since row phases are adjustable, this is equivalent to the TBM 
in (2). 

(d) In conclusion, as long as F e C4, the full-mixing matrix is equivalent to the TBM 
in (2), characterized by the physical parameters (;^, :^,0). 



2. 



F e C5 . Any F in this class obeys = 1, hence its possible eigenvalues are 
1,+i, — i,— 1. Since we want it to be non- degenerate, has determinant +1, and 
since its character according to Table 1 is real, we conclude that its allowed eigen- 
values are each occurring once and only once. For definiteness, we will let 
F — aib — (1342) := Ci. This is completely general because the four permutation 
objects in 5*4 can be assigned any label we wish to. 

Once again, G must be chosen so that {F, G} generate S4. To see how to find such a 
G, write the presentation of 5*4 in (4) in a slightly different form, 5*4 — {a, c|a^ = — 
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(ac)^ = 1), which shows that the generator a is of order 2, the generator c of of order 
4, and ac is of order 3. Trying all 9 possibilities in C2 and C3 with c — ci — (1342), we 
see that only the following four a's, ai = (12), a'^ = (34), a'^ = (13), and 03 = (24), can 
yield an order-3 ac. Thus G is confined to one of these four cases. As a consequence, 
the allowed triplets are G = Ai and ^43. 

In order to find out what these matrices are in the F-diagonal basis, we first write 
them in the 6-diagonal basis, then convert them by a similarity transformation to the 
F-diagonal bases. 

In the 6-diagonal basis, 

/ 1 
-2 



-2a; -2a;2 \ 



V-2 



-2uj 



This F can be diagonalized by the unitary matrix 

/ UJ —UJ 
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to yield V^FV = diag(l, i, -i) := F' . 
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(a) G = Ai . In the F-diagonal representation, [Gi, 6*2! ^3] 
y^[Gio, G30, G2o]V. The full-mixing matrix is then given by 
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where Uq is the TBM mixing matrix given in (2) with the second and third 
columns switched. We do that switching to put U into a more familiar form with 
a vanishing reactor angle. 

Multiplying the first, second, and third columns by g'^Ve^ g-7r4/i2^ ^5^/12^ ^^^q 
first and second columns by i, —1, we can turn U into the equivalent form 

/1/V2 1/V2 \ 

-1/2 1/2 I/V2 , (9) 
V 1/2 -1/2 I/V2) 
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whose mixing angles are given by the physical triplet ^i^) ■ In this case, the 
reactor angle remains to be zero, but both the solar and the atmospheric mixings 
are maximal. 

(b) G = ^3 - First, note that = [(24), (13), (24) (13)] = Ci[(34), (12), (34)(12)]cr\ 
where ci = (1342). In the 6-diagonal representation, [03,03,03] = 
ci[ai,a\,a'l]cl^ = F[Giq,Gzq,G2q]F~'^ . In the F-diagonal basis, [^1,(^2,^3] 
[G',G,G"] = V^F[Gio,G^o,G2^]F-^V = F'V'^[G^^,G^^,G2o]VF'-\ The full- 
mixing matrix is therefore F'V^U'rj.^^ = F'U^, with f/5 being the U in (8). Since 
F' is a diagonal phase matrix, the full-mixing matrix here is equivalent to (8) 
and (9). 

(c) In conclusion, as long as F e C5, the full- mixing matrix is equivalent to U of (9), 
characterized by the physical parameters i:^^ 

To summarize, up to equivalence, the full-mixing matrix of 5*4 is either given by the TBM 
in (2), characterized by the physical parameters (^, 0), or the matrix (9), characterized 

by (75' 73' 0)- 
B. A4 

The character table of A4 is given in Table 2. 





Ci 


C2 


C3 C4 




(1) (12)(34) (123) (132) 


\C^\ 


1 


3 


4 4 


x' 


1 


1 


1 1 


x' 


V 


1 




x' 


I" 


1 




x' 


3 


-1 












h 



Table 2. Character table of A4 



In this case the order-3 elements are divided into two classes, C3 and C4, with the elements 
in C4 being the square, or the inverse, of the elements in C3. Since the generator b is of order 
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3, it is either in C3 or C4, which differ only by the arbitrary naming of the four objects. For 
definiteness we will choose it to be in C3. Since its three dimensional character vanishes, 
up to permutation of the entries, the diagonal form of b is again b — diag(l, cu, a;^). For 
definiteness, once again we will let b = (134). 

The generator a is of order 2, so it must be in €2- To be a generator, it must contain the 
permutation of the object '2', but then all of them do. Hence there are three possibilities: 
ai = (12) (34), 02 = (13) (24), 03 = (14) (23). These three mutually commute, and the 
product of any two is the third, so they can be taken to be G, G", G" respectively. 

In the diagonal form, a non-degenerate F must he F — b — diag(l, a;, a;^) — Fq, up to 
permutation of its entries. We must now find out the representation of A = (01,02.03) in 
this basis. 

We know from (3) and the discussions below that equation that Fq and G20 generate A4, 
hence we may take 02 = G'2o- Then oi = (12)(34) = (134)(13)(24)(143) = ba^b^ := G[, 
and 03 = (14) (23) = b'^aib :— G3. Since the invariant eigenvector of G2 is proportional to 
(1, 1, 1)^, and invariant eigenvectors of G[ and G'^ are respectively Fo(l, 1, 1)^ = (l,a;,a;^)^ 
and Fq^(1, 1, 1)^ = (1, cu^, cu)^. Consequently, the full- mixing matrix of A4 is the Cabibbo- 
Wolfenstein mixing matrix 



/I 1 1 \ 



1 



OJ 1 oj"^ 



(10) 



u 



V3 



Vu;^ 1 00/ 



whose triplet of physical parameters can be shown to be ( 



1 1 e'""/^ 
a/2' V2' a/3 



) 
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C. As 



The character table of A5 is given in Table 3, 







C2 


C3 


C4 


Cs 




(1) 


(12)(34) (123) (13524) (12345) 




1 


15 


20 


12 


12 


x' 


1 


1 


1 


1 


1 


x' 


4 





1 


-1 


-1 


x^ 


5 


1 


-1 








x' 


3' 


-1 





-l/ip 


V 


x^ 


3 


-1 













a 


h 


ab 


{abf 



Table 3. Character table of 



where ip is the golden ratio (1 + V5)/2 = -2cos(47r/5) = 1.618, and l/^p = (-1 + V5)/2 = 
if-l = 2cos(27r/5) = 0.618. 

The elements in class C5 are the square of the elements in C4. From the presentation of the 
group given in (4) , we see that the generators a and b must in classes C2 and C3 respectively, 
and ab must in C4 or C5. Which of the two it belongs simply depends on how we label the 
five objects in ^45. 

For definiteness, let b = (134) as before. To be a generator of ^5, a must contain the 
numbers '2' and '5' in separate cycles, such as 013 = (12)(35) or = (15)(23). These 
two commute and produce a product a"^ = (13) (25), which incidentally is not a generator 
because '2' and '5' are in the same cycle. In addition to Ais — [ais,a[^,a'l^], there are 
two other commuting triplets, A14 — [ou, a'14, a'/J and ^434 = [034,034,034], with 014 = 
(12)(45),a;4 = (15)(24),<4 = (14)(25), and 034 = (32)(45),a^4 = (35)(24), a^'4 = (34)(25). 
In each case, a or a' is a possible generator of ^Is, but not a". Any of these three triplets is 
a possible choice for G = [G, G', G"]. 

F may be in C3 or C4, but we will only consider the case F G C3 here [12] . As usual, 
without loss of generality we may choose F — b — (134) = diag(l, o", o"^) = Fq as before. 
The question is what is the three-dimensional representation for the generator a. 

Representation of can be found in [5], but unfortunately it is not in the 6-diagonal 
basis that we need. Nevertheless, from that representation, one can work out the defining 
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representation a in the 6-diagonal basis to be 



V5 
6 



^{-2V5 + 2iV3- iVbVS) 



12 



1 VE 

2 6 



One can check that — — 1 and (a6)^ = 1, so these a and 6 are indeed generators of A^. 

This a can be equated to any one of the generators in C2. For definiteness, let it be 

= = (12) (35). We must also find the matrix representation of its commuting partners 

^'13 •= '^13 •= They are in the same class as a, so each of them can be written in 

the form gag~^, for some g which itself a product of the generators a and b. It can be shown 

that g = {bafb'^ = (125) and g'^ = b{baf = (152) convert a = (12) (35) into a" = (13) (25). 

Using the representation of b and a, one can calculate from (11) to obtain 

/-8 \ 

a" = l 1 + 3^5 + V3i(l-V5) , (12) 

8 

V 1 + 3V5 - V3i{l - VZ) / 

and then from a' — aa" that 



1 _i_ vl 

2 6 



^- J^(^/5-^V3-^V5V3) 
V3- n(\/5 + iV3 + i^y3) i + ^(-2y5-2iV3-iV5V3) 



^(-2V5 - 2iV3 + iV5V3) 



V5 
6 



(13) 



The unitary matrix to diagonalize simultaneously a, a', a" so that WaU 
diag(+l,-l,-l), Wa'U = disig{-l, +1,-1), [/V?7 = diag(-l, -1, +1) is 

( ^ f-^ \ 



73 



i(l+a;(^)(^-i Kl+w^)^ ^i(l + a;(^) . (14) 
V|(l+wV)v-' -Kl+wV)^' ^i(l + a;V)/ 
Writing (1 +ujLp) = (l + cj^y?)* = \f2e^^, then multiply the second and the third rows of U by 
e~*^ and — e*^ respectively, and the third column by —i, we can turn U into the equivalent 
form 



U 



1 



V2ip 
1 

V2ip 





\ 




73 




V2 




\/3 , 
V2/ 


V2 



(15) 
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whose physical parameters are given by the triplet (;^^);^)0). This is the full-mixing 
matrix of A5 if (G, G', G") = (013, a^g, a'/g) = (a, a', a") = ((12)(35), (15)(23), (13)(25)). 

Next, consider the other triplet (G, G", G") = (om, a[^, a'/4) = ((12)(45), (15)(24), (14)(25)). 
Let g = {baf = ((134)(12)(35))3 = (15243). Then {au, a[^, a'l^) = g{a, a" , a')g-\ Note that 
the positions of a" and a' are reversed. The full-mixing matrix for this triplet choice is 
therefore U — XU, where U is the matrix in (14) with the second and third columns 
interchanged, and 

/ 4 1 - iVS + V5{3 + iVS) l + iV3 + V5{3-iy/3)\ 

4(1 -iV3) 1 + iVS - V5{3 - iVS) l-iV^VS . (16) 

V4(l+iV3) l + iVsVS 1 - i\/3 - V5(3 + iV3) y 



X = (baf = — 
^ ' 12 



The result is 



V3 



1 



— GOP"*'' J-co-ip-*'' —jy^p-iv 



(17) 



where 77 is defined by ^ =1=1 = 75 

Multiply the second row by e*^, the third row by — e~"', the second column by —1, the 
third column by i, and interchange the first and second columns, we get back to the mixing 
in (15). 

Finally, similar calculation shows that if (G, G', G") = (034, 034, 034), then its full-mixing 
matrix is also equivalent to (14). 

In summary, as long as F G G3, the full-mixing matrix of is equivalent to (15), a 
mixing which is characterized by the triplet physical mixing parameters {-^^1 0)- 

IV. CONNECTION BETWEEN GROUP THEORY AND DYNAMICS 

Charged-lepton masses come from Yukawa interactions in the Standard Model. If neu- 
trino masses arise from the type-1 seesaw mechanism, then the Lagrangian responsible for 
flavor structure of mass matrices can be symbolically written as 

L = -c^eReL(f - h'^NRPLX^ - ^H^N^Nrx^ + h.c, (18) 

where e^, cr are the left-handed and right-handed charged leptons, v^, Nr are the active 
and heavy Majorana neutrinos, and (f)^ tX'^ ^-re the fiavons. Isotopic spin and spacetime 
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details are omitted, so this Lagrangian should only be used to discuss the flavor structure 
of mass matrices. In that regard every term is assumed to be invariant under a horizontal 
symmetry group Q, with c^,h^, h'^ being the Yukawa coupling constants. The index 7 labels 
the irreducible representations (IR) of Q, and a sum over all 7 is understood. The IR of the 
fermions are not explicitly specified. 

For type-II seesaw, we replace the last two terms with a single term, hence 

C = -c^e.ReL(tP - K^l^lX^ + h.c. (19) 

To obtain the mass matrices from (18) or (19), vacuum expectation values of the flavons 
are introduced to break the ^-symmetry. In order for the residual symmetries to be pre- 
served after the breaking, F in the charged-lepton sector and Gi in the neutrino sector, the 
vacuum alignments have to be invariant eigenvectors (eigenvector with eigenvalue +1) of 
the corresponding residual operators [3], 

fW(</,7) ^ (<^7)^ Gt'\x')-{X'), Gt'\x')-{X'), (20) 

where F^^^ and G\"'^ arc the representations of F, Gi in IR7. If we want to obtain a full- 
mixing matrix, then (20) is required to be obeyed for i = 1, 2, 3, but if we want to construct 
a partial-mixing matrix, then the equation has to be satisfied only for that i. 

Assuming the invariant eigenvectors to be unique, it follows from (20) that (x''') = ix'^) 
up to normalization, so we will set x'"' = X^ in (18) from now on. 

Note that (20) does not rely on the specific form of (18) or (19). All that it requires is 
the neutrino mass terms to be distinct from the charged-lepton mass terms, and that each 
has its own valon. (20) would be equally valid if valon fields appear quadratically in the 
Lagrangian, for example. 

Our main goal is to keep the effective left-handed mass matrices and M^j invariant 
under the residual symmetry. This can be accomplished by keeping the broken Lagrangian 
invariant, as is done above, but in the case of type-I seesaw, since the Dirac mass matrix 
My enters quadratically in My = M^ M]^^My, where Mm is the Majorana mass matrix of 
Nr, invariance of My can also be preserved indirectly [6] if the last term in (20) is replaced 

^yGf\x')^-{x'). 

Once we know the IR of the residual symmetry operators, vacuum alignments can be 
calculated. With those alignments, one can construct many models that automatically give 
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rise to the desired neutrino mixing. Yukawa coupling constants are used to fit the fermion 
masses, and in the case of partial mixing, also the remaining mixing parameters. Different 
models are obtained by choosing different IR assignments of e^, cr, vl, Nr, 0, and x, though 
for a successful model we must make sure that enough terms are present in each sector to fit 
all the masses and the remaining mixing parameters, in spite of the zero vacuum alignments 
forced upon us in those IR without an invariant eigenvector. 

Let us see how this works explicitly. Let {aa,/3b\^c) be the Clebsch-Gordan (CG) coeffi- 
cient of the group Q, coupling orthonormal states \aa) in IRa with the orthonormal states 
\/3b) in IR^ to get an orthonormal state |7c) in IR7, 

a,b 

\aa)\(3b) = Y.\jc){jc\aa,(3b). (21) 

7,c 

The CG coefficients must obey the group-invariant condition 

c' a',b' 

^gi'${jc'\aa,pb) = T^Maa' , pb')g^lgi^l (22) 

c' a',b' 

for every e ^, or at least for every generator of ^. It is important to note that if we change 
the IR bases of g, the CG coefficients may have to be changed as well. In the case of S'4, the 
CG coefficients to be used will be those that correspond to the IR given in Table 4. 

Suppose eR and belong to IRa and lRf3 respectively, and the fiavons belong to IR7. 
Then the charged-lepton mass matrix can be read out from either (18) or (19) to be 

(Me)„, = Ec7(«*«,/3&|7c)(</'^). (23) 

The invariance of the charged-lepton mass matrices Me and Me :— M^Me under the residual 
symmetry F can be seen from (20) and (22): 

i7^(")t7^^^^(/3) ^ F^^^^MeF^^^ = Me. (24) 

Similarly, if Nr, i/l belong to IRa and lRf5 respectively, then the neutrino mass matrices 

are 

{M,)ab = Y.h'-y{aa,(3b\jc){x2), {Mr,)aa' = ^ ^7(«a> «a'l7c)5(x7), = MJMj^1M,(25) 

7,c 7,c 
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where (aa, aa' 170)5 is the CG coefficient symmetric in a and a'. Once again the invariance 
relations 

follow. Note that because of the Majorana nature of the neutrino, we have assumed G^-"'^^ 
to be real. For indirect models, the first equation in (26) should have a minus sign. 

In the rest of this section, we will illustrate these operations by considering the direct 
full-mixing models of the groups Q = A4 and S4. In the case of A4, partial- mixing models 

arc also discussed because they can give give rise to the phenomenologically interesting 
tribimaximal mixing. 

The first task in each case is to find out the IR of the residual symmetry operators. 



Remember from Sec. IV that there are two full-mixing modes, one giving rise to the TBM, 
and the other yields a mixing with too large a solar angle. We shall consider how dynamical 
models can be constructed in each case, starting with those giving rise to the TBM. 



1. 



F — Fq, [Gi, G2, G3] — [Gio, G20, G; 



Table 4 shows the IR and invariant eigenvectors of the residual symmetry operators 
in the basis where F is diagonal. 



IR 


1 


1' 


2 


3' 


3 


F 


1 


1 


diag(a;, cu^) 


diag(l, u, u"^) 


diag(l, cu, UJ^) 


(0) 


1 


1 


(0,0) 


(1,0,0) 


(1,0,0) 


Gi 


1 


-1 


cri 


— Gio 


Gio 


G2 


1 


1 


diag(l, 1) 


G20 


G20 


G3 


1 


-1 




— G30 


G30 


(x) 


1 





(1,1) 


(1,1,1) 


(0,0,0) 



Table 4. Irreducible representations and invariant eigenvectors of the first set of S4, 

residual symmetry operators 
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The one-dimensional IR 1, 1' in Table 4 are obtained directly from Table 1 because in 
that case character is the IR. The three-dimensional representations 3, 3' are copied 
from those given in the previous section. For 2, since — 1 and its character from 

Table 1 is Tr(C4) = —1, its diagonal form must be diag{u, u"^). The other generator Gi 
must not be diagonal, or else the 2-dimensional representation is not irreducible. Since 
Gf — 1, within equivalence we may choose it to be the Pauli matrix c''! = ^ ^ q^- 

With these IR of the residual symmetry operators, it is straight forward to work out 
the invariant eigenvector (0) of F, and the simultaneous invariant eigenvector (x) of 
Gi,G2, and G^. The column vectors are displayed in the table as row vectors for 
printing convenience. Some of the vacuum alignments are zero either because does 
not have +1 eigenvalue, or because there is no common eigenvector of G", G2 , G2 all 
with eigenvalue +1. 

With the help of Tables 1 and 4, many direct models possessing TBM can be con- 
structed. These models differ from one another in their IR assignments for the right- 
handed fermions and flavons. 

We saw in the last section that the residual symmetry operators for ^'4 and A4 are 
in IR3, hence the left-handed fermions must transform as 3. One can then work out 
from the character table the allowed IR of the fiavons for every assignment of the 
right-handed fermion. The result is listed in row 2 of Table 5. Since (0) = for IR2 
and (x) = for IR3, neither of them can be used to construct models. Taking that 
into account, the final number of free parameters available in model construction is 
listed in the third row for the charged-lepton sector, the fourth row for the neutrino 
sector if it is given by a type-II seesaw, and the fifth row for the neutrino sector given 
by a type-I seesaw. Those in the neutrino sector require a little bit of explanation. 

For type-II seesaw, ul is coupled to itself so only 3x3 couphng is relevant. For 
type-I seesaw, the heavy Majorana neutrino Nji can be assigned to any IR. The first 
of the two numbers in the last row refers to the number of free parameters in the 
Dirac neutrino mass term. Once the assignment a of N^i is determined, the number of 
parameters in the Majorana mass terms is determined hj a x a, and that appears as 
the second number in the last row. If we consider indirect models [6] , then the number 
of parameters is even larger than that. 
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RHF 


1 


v 


2 


3' 


3 


flavoiis 


3 


3' 


3' . 3 


1'. 2.3'.3 


1.2 . 3' . 3 


e param 


1 


1 


2 


3 


3 


typell v 










3 


typel V 


0+1 


1 + 1 


1 + 3 


3 + 3 


3 + 3 



Table 5. Number of free parameters in 5*4 models with a TBM mixing 

Let us see what Table 5 is telling us. First look at the charged-lepton sector, in row 
3. Since the number free parameters always happens to match the dimension of the 
representation, there are always just enough of them to fit the three charged-lepton 
masses. 

For type-II seesaw, there is again exactly three free parameters to fit the three active 
neutrino masses. For type-I seesaw, the number of free parameters shown is the sum 

of the My and parameters, which is generally more than enough to fit the active 
neutrino masses. The remaining ones can be used, for example, to tune the properties 
of the heavy neutrinos. 

As a concrete illustration, we will write down the mass matrices of a type-I model 
where e/,, cr, vl, all belong to 3. It is the kind of assignment that a 5'O(10) grand- 
unified theory would want. In that case, according to Tables 4 and 5, Mg, M^, and Mat 
should contain exactly three parameters each. Using the CG coefficients in the Fq- 
diagonal basis, these mass matrices can be worked out to be Mg = diag(ci + 2C3/, ci — 
C3' + C3, Ci — C3/ — C3), where Ci = Ci/-\/3, c^i — C3//3-\/2, and C3 = Cs/y/E. Furthermore, 





/hi + 2/13/ 




- hy 


h2 




Mn = 


h2 - 


h2 


■1 r 


hi 





(27) 



\ /l2 - hy hi + hy h2 - hy J 



where hi — hi/y/3, hy — /i3'/3-\/2, and h2 — h2/VE- The expression for M^, is 
identical except that h^ is to be replaced by h'^. 

This example illustrates some general points discussed above: 

(a) Each of Mg, M^,, is specified by three independent parameters, as shown in 
Table 5 above. 
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(b) The three parameters in Mg can be fitted by the three charged-lepton masses. 

(c) There are six parameters to specify Mj, = MjM^^Mj,. In additional to fitting 
the three active neutrino masses, we can use the remaining three to specify the 
heavy Nr masses. 

(d) Me is automatically diagonal. This is because we are using CG coefficients in the 
Fo-diagonal (hence Mg-diagonal) basis. As remarked below (22), if we use some 
other basis, then the CG coefficient will be rotated and the resulting may no 
longer be diagonal. 

(e) Whatever the parameters are, both M,, and turn out to be magic (all 
row sums are equal) and 2-3 symmetric (invariant under a simultaneous ex- 
change of the second and third columns, and second and third rows). Hence 
M^, = M^^ Ml, is also magic and 2-3 symmetric, and the mixing matrix U is 
automatically TBM [7]. 



2. 



F — GiqFq, G2, G3] — [Gw, G20, G 



f30 



This is the other allowed full mixing of S4, given in (6), (7), (8), and (9). Unlike 
the TBM mixing, whose models abound in the literature, this one is seldom discussed 
presumably because its phenomenology is less attractive. Another reason might be 
that the appropriate vacuum alignments in this case are a bit strange, as we shall see. 
Nevertheless, we can follow the general procedure to construct models that lead to 
this kind of mixing. 

There are two different ways to proceed. Either we go directly to the F-diagonal 
basis, as prescribed before, or we stay in the Fo-diagonal basis that we are familiar 
with. The former has the advantage of having Mg automatically diagonal, but the 
CG coefficients in the F-diagonal basis must be worked out anew. If we travel the 
latter route, the CG coefficients are the same ones used in (27), but we have to work 
out the new alignments, and have to perform the diagonalization of Mg afterward. 
To illustrate better the connection between the group-theoretical and the dynamical 
approaches, we choose to work with the latter. 

In the 6-diagonal basis, the IR of Gj is the same as in Table 4, and the IR of F can be 
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obtained from Table 4 by calculating GiF. The resulting IR as well as the invariant 
eigenvectors are shown in Table 6 below. Since the IR for Gi remains unchanged, that 
part of the table is identical to Table 4. The only difference comes in the upper part 
of the table in which F is replaced by GiF in every IR. 



IR 


1 


1' 


2 


3' 


3 


F 


1 


-1 


(Tidiag(a;, a;^) 


—GioFq 


GiqFq 


(0) 


1 





(1,0;) 


(0,0,0) 




Gi 


1 


-1 


o-i 


—Gio 


Gio 


G2 


1 


1 


diag(l,l) 


G20 


G20 




1 


-1 




—Gso 


G30 


ix) 


1 





(1,1) 


(1,1,1) 


(0,0,0) 



Table 6. Irreducible representations and invariant eigenvectors of the second set of S4 

residual symmetry operators 

The parameter count for this kind of models is shown in Table 7. The only variation 
from Table 5 comes from the different vanishing (0)'s, now occurring in 1' and 3' 
instead of 2. 



RHF 


1 


1' 


2 


3' 


3 


flavons 


3 


3' 


3', 3 


l',2,3',3 


1,2, 3', 3 


e param 


1 





1 


2 


3 


typell v 










3 


typel z/ 


+ 1 


1 + 1 


1 + 3 


3 + 3 


3 + 3 



Table 7. Number of free parameters for the other 5*4 model 

As a result, in order to have enough parameters to fit the charged-lepton masses, cr 
must now be in 3, not 3', nor a combination of 1 and 2. Other than that, things are 
pretty much the same as before except that is now no longer diagonal. To see the 
difference between this case and TBM, let us illustrate again with the model in which 
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eL,eR,i^L, Nr all belong to 3. In that case, using the Fo-diagonal CG coefficient as 
before, but now a different (0), we get 



(28) 



C2 + co'^cs/ a/3 Ci - UJC3/ ^/3 UC2 

\ UJC2 - C-il "\/3 C2 Ci + UJC-il a/3 / 

where c\ = ci/^/S, C2 = C2/\/6, and C5 = c^l^2. Since the (x)'s are the same as in 
Table 4, the neutrino mass matrices My and Mat are identical to those given in (27), 
which are diagonalized by the TBM matrix Uq of (2). If for convenience we switch the 
second and the third columns, then it is C/q as explained below (8). This U'^ however 
is not the mixing matrix because Mg is not diagonal. 

Since Mg commutes with F, which is non-degenerate, the matrix V in (7) that diago- 
nalized F must also diagonalize Mg. Indeed, 

/ci + 2a;2c2 \ 



Ci - U?C2 + WC5 



V 



(29) 



ci - a;"C2 - a;c5 / 

which gives rise to the charged-lepton mass-squard matrix V^M^V — 
diag(mg, m^, m^), with 

ml = {ci - €2)^ + 3cl, ml.^ = cl + cl + cl + CiC2Tici-C2)c5. (30) 

The neutrino mixing matrix is now given hy U = V^^Uq, whose explicit result is already 
given in (8) and (9). 



B. A4 

We shall discuss for this group both the full-mixing and the partial-mixing models. For 
full- mixing, (x) is the simultaneous invariant eigenvector of Gi, 6*2, and G3, and is denoted by 

iXaii)- For partial mixing defined by Gi, (x) is the invariant eigenvector of that Gi, denoted 
as (xi). Table 8 exhibits the IR of the residual operators and the vacuum alignments. 
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F0G20F0 := 6*30 


(X3) 
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1 


1 


(1,^2,^) 


iXall) 


1 


1 


1 


(0,0,0) 



Table 8. Irreducible representations and invariant eigenvectors of A4 
The parameter count for these models is given in Table 9. 
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3 + 3 



Table 9. Number of free parameters for the full- mixing and partial- mixing A4 model 

We see once again that no matter how cr is assigned, there are always just enough 
parameters to fit the charged-lepton masses. For type-II neutrinos, after using three to 
fit the neutrino masses, there is one left over to fit the remaining mixing. It is this extra 
parameter that can turn a partial-mixing mode of into a TBM, thereby promoting the 
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A4 symmetry to a S4 symmetry. In the case of full mixing, listed under 'all', the mixing 
matrix is completely determined to be the Cabibbo-Wolfenstein matrix, so the number of 
parameters is just enough to fit the neutrino masses. For type-I mixing, as before, there are 

more parameters than necessary to fit the active neutrino masses, and the remaining mixing 
if in the partial- mixing mode, so the extra ones can be applied to tune the properties in the 
heavy neutrino sector. 

V. HIGGS POTENTIAL OF A4 

We established in Sees. II and III a connection between the mixing matrix and horizontal 
symmetry, purely on symmetry grounds, without the presence of Higgs nor a Lagrangian. 
Wc also saw in Sec. IV that if fiavons are introduced in a Lagrangian, then their vacuum 
alignments must be invariant eigenvectors of the residual symmetry operators for the two 
formalisms to agree. In a dynamical model, vacuum alignments are obtained by minimizing 
a Higgs potential, so for consistency we must find a suitable Higgs potential to yield just 
those alignments. 

A completely generic Higgs potential contains so many ^-invariant terms that just about 
any vacuum alignment can be obtained by adjusting the self-coupling parameters. In order 
to confine the solutions to invariant eigenvectors, a rationale must be found to chop down 
the number of terms in the potential. The usual way to accomplish that is to impose an 
additional symmetry on top of Q. 

Take Q — A4 for example. It contains four irreducible representations, 1, 1', 1", 3. Since 
the only non-trivial alignment occurs in 3, let us consider a renormalizable Higgs potential 
V{iIj, ip*) of the complex triplet flavon field ip = (oi, 02, 0.3)"^. It turns out that if V contains 
only and ip* ip* ipip terms, without terms like ipip, ip* ipip , 'ip*ipipip. etc., then its minimiza- 
tion would yield just the invariant eigenvectors of A4. Such a potential can be arrived at 
by assigning ip a non-trivial U{1) quantum number, ip* the opposite quantum number, and 
imposing an additional U{1) symmetry on V. A suitable Z{N) would do as well. 

In addition to showing that the solution to its equations of motion are invariant eigen- 
vectors, we shall also show that alignments for A4 members of the same class have the same 
Higgs energy. These general properties remain valid if A4 is replaced by any simply reducible 
group Q, as we will see in the next section. 
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Such a potential has been considered before [8] . For the quartic terms, if we choose to 
couple the two ■^'s first, and separately the two then we obtain four terms depending 
on whether V'V' are coupled to 1, 1', 1", or 3. Using V"^ to denote the projection operator 
for irreducible representation a, the most general V is 

V = {^hV^ + fvV + fv^V" + hV^\^^) - IJ?{^\^), (31) 

where fa are the coupling constants. 

With the help of explicit A4 CG coefficients given in Appendix A, V can be written as 

+ /I'bif + /i-bi'f + 2/3(A;' - \yz?)) - 11% (32) 

where 

yi = + 20203, yv = 03 + 2aia2, yv ^ al + 20103, 

yz = 0102 + 0203 + 030^, A; = |oi|^ + |o2|^ + |o3|^. (33) 



The equations of motion are obtained by setting dV/ dai = and dV/ da* = 0. In this way, 
we obtain the following equation and its complex conjugate: 
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{^3) 













- r 



V«3/ 



(34) 



It is much easier to solve this cubic equation of Oj by first converting it to a quadratic 
equation in ya and k. To accomplish that, we multiply (34) on the left respectively by 
(0^,0^,0^), (0^,0^0^), (0^,0^, a:^) to get 



/ilyir + /I'lyi'l' + /i"l?/i'f + 2/3(A;' - |z/3r) - = O, 

fiywl" + fi'Vi'Vl + fvyvyv + f^ikyl - yl) - (3/2)^^1/3 = 0, 

fmyl' + fvyvyl" + fvyvyl + f^iky^ - yf) - (3/2)//^y3 = o. 

Substituting (35) into (32), we can simplify the Higgs energy expression to 

V = -\^?k. 



(35) 
(36) 
(37) 

(38) 
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Subtracting the complex conjugate of (37) from (36), and ask the equahty to be true for all 
fa, we obtain 



ywl" ^ vlvv ^ yvVi" ■ (39) 

Let us now proceed to solve (39). 

1. If yi = 0, then yi' or yi'/ must vanish. More generally, if one of the ys (s = 1, 1', 1") 
vanishes, then one or both of the other two must also vanish. In this way, we obtain 
four possible solutions: y :— [yi, yy, yi"] oc [0, 0, 0], [1, 0, 0], [0, 1,0], [0, 0, 1]. 

2. If none of the ys vanish, then it follows from (39) that \yi\ = \yii\ = 
Dividing (39) by yiy^, we get {yv/yi)* = {yi'/Vi) = {Vv I yi){yv / Vi)* ■ Hence 
{yvlVif = 1, {yvlvif = 1, and (yv/yi) = {yvlvif- The solutions are therefore 
y = 1, 1], yi[^,^,^\ and 

Substituting these solutions of yg into (33), we can obtain the solutions of = (oi, 02, 03). 
Because of the U{1) invariance, the solutions below can always be multiplied by an arbitrary 
phase factor of unit modulus. 

1. (a) Suppose yv — yyi — 0. 

i. One obvious solution is a2 = 03 = 0. In that case, ai = if and only if 

yi = o. 

The solution = (0, 0, 0) gives rise to A; = ^3 = 0, and a Higgs energy V — Q 
from (38). Obviously it also satisfies the equation of motion (34). 

The solution if) = ^/k:{l, 0, 0) gives rise to 7/3 = and yi = k. It satisfies (34) 
and (35) if A; = 3/x^/2(/i + 2/3). The Higgs energy is then given by (38) to 
bey = -3//72(/i + 2/3). 

ii. There is also a solution with 02 7^ 0. In that case, yi' — yi" — implies 
ai = — a|/2a2 and = a|/a2. The latter yields 03 = a2Uj"^ for m = 0, 1, or 
2, and the former demands ai = —a2(^^"''/2. From (33) we now obtain yi = 
9a|cc;™/4, k = 9|a2p/4, and y^ = 0. We can now check that (34) is satisfied 
provided k = 3/^72 (/i + 2/3), so the Higgs energy V = -3/^72 (/i + 2/3) is 
degenerate with that of the i/j — Vk{l, 0, 0) solution. 
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This solution (actually three solutions corresponding to m — 0,1,2) can 
then be written as t/j — (-\/^/3)(l, — 2a;"*, — 2a;^"*), or more explicitly, i/j — 

(v^/3)(l, -2, -2), (v^/3)(l, -2a;, -2u^), (v^/3)(l, -2a;2, -2a;). 

(b) The solutions for yi = yyi = 0, yy ^ can be similarly obtained. It leads 
to the solutions i\) = ^(0, 1, 0), {.Jk/?,){-2, 1, -2), (V^/3)(-2a;, 1, -2a;2), and 
(-\/^/3)(— 2a;^, 1, —2a;). The values for k and V are respectively k — 3/i^/2(/i/ + 
2/3) andy = -3Ai72(/i' + 2/3) 

(c) The solutions for yi = yy = 0, yy, are ip = y/k{0,0, 1), {^fk/3){-2, -2, 1), 
(V^/3)(-2a;,-2a;2,l), (V^/3)(-2a;2, -2a;, 1). The values for k and V are re- 
spectively k = 3//V2(/i» + 2/3) and V = -3//72(/i" + 2/3). 

(d) In short, the solutions of (b) and (c) are obtained from (a) with appropriate 
permutations. 

2. (a) y = yi[l, 1, 1]. From (33), we need to have + 20203 = 03 + 2aia2 = a| + 20103. 

The first equality leads to (ai — 03) (ai + 03 — 202) = 0, and the second equality 
leads to (02 — 03) (02 + 03 — 2ai) = 0. The solution is ai = 02 = 03, hence 
-0 = {^^k/3){l, 1, 1). Using (33), we get yi = = yi" ^y^^k. 

We can check that (34) is satisfied provided k — 3/x^/2(/i + /i/ + fyi). The Higgs 
energy is then V = -3iJ,^/4:{fi + fr + /i//). 

(b) y — yi[l,a;,a;^]. Prom (33), we need to have al + 20203 = <^^(ai + 201O2) = 
a;(o| + 201O3). The first equality leads to (oi — a;o3)(oi -|-a;o3 — 2a;^02) = 0, and 
the second equality leads to (02 — a;^03)(o2 + a;^03 — 2a;oi) = 0. The solution is 

Oi = a;^02 = a;a3, hence ip = (-^/c/3)(l,a;,a;^). Using (33), we get yi = oj'^yv = 
<^yi" = y3 = k. 

We can check that (34) is satisfied provided k — 3/i^/2(/i + /i/ -I-/1"). The Higgs 
energy is then V = -3//^/4(/i + fy + /i//). 

(c) y = yi[l,u'^,u]. The solution is = {^k/3){l, u"^, u), with k = 3/x^/2(/i + fii + 
/rO- The Higgs energy is t = Si^V^f^ + f,, + f,,,). 

3. In summary, the solutions can be classified according to their energy into five cate- 
gories: 
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(a) y = 0. In this case ^ = (0, 0, 0). 

(b) V = -3//V4(/i + 2/3). In this case V = v^(l,0,0), (\/fc/3)(l, -2, -2), 
(\/^/3)(l, -2a;, -2a;2), or (\/fc/3)(l, -2a;2, -2a;), with k = 2>iJ?j2{fi + /a). 

(c) V = -3//74(/i, + 2/3). In this case -0 = v^(0,l,0), (v^/3)(-2, 1, -2), 
(v^/3)(-2a;, 1, -20;^), or (^/^/3)(-2a;^ 1, -2a;), with k = 3//V2(/r + /a). 

(d) \/ = -3/i74(/i» + 2/3). In this case ^ = v^(0,0, 1), (v^/3)(-2, -2, 1), 
(V^/3)(-2a;, -2a;2, 1), or (V^/3)(-2a;2, -2a;, 1), with k = "^ii^ I2{jv' + /a)- 

(e) y = -3//V4(/i + jv + /i")- In this case i) = (\/^/3)(l, 1, 1), (v^/3)(l,a;,a;2), 
or (v^/3)(l,a;2,a;), with k = 3/xV2(/i + fr + /i")- 

In particular, if fi, fv, fv, fs are all positive, energetics dictates that the solution ip — 
(0, 0, 0) is disfavored. This is the solution that leads to the A4 full-mixing matrix of Cabibbo 
and Wolfenstein, rather than the experimentally correct mixing of the tribimaximal type, 
obtainable using the first solution of (b) in the charged lepton sector, and the first solution 
of (e) in the neutrino sector. 

Moreover, if fi < fit and /i", then it is (1,0,0) etc. that is energetically favored, and not 
(0,1,0), (0,0,1), etc. 

In this way, we see that there is a rational why A4 might prefer partial mixing, which may 
lead to a tribimaximal matrix, rather than a full mixing that leads to a Cabibbo- Wolfenstein 

matrix. 

We proceed now to discuss the connection between these vacuum alignment solutions and 
the residual symmetry operators. Table 10 shows the vacuum alignment solutions above, 
and the operator for which the the alignment is an invariant eigenvector. F, G'^, 6*2, G3 are 
the generators of ^4 in Table 8, but with the subscript omitted. 
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alignment 


operator 


(1,0,0) 

(1,-2,-2) 
(1, -2cj, -2a;2) 
(1, -2w^ -2uj) 


FG2F, F^GsF^ 

F^G2, G2F 

G2F'^ , FG2 


(1,1,1) 
(1,0;, ^2-) 

_L . uJ . uJ ^ 


G2 

G[ = F^G2F 

G', = FG2F'- 


(0,0,0) 


{G'l, G2, G3) 



Table 10. Dynamical alignments are invariant eigenvectors of appropriate operators 



Here are some remarks concerning Table 10. 

1. The first seven solutions in the left column are the invariant eigenvectors of the 11 
operators in the right column, each of them responsible for a partial-mixing. Together 
with the identity operator, these twelve operators together form the A4 group. 

2. The last solution is the simultaneous eigenvector of G[, 6*2, and G'^ relevant for full- 
mixing. 

3. Solutions under 3(c) and 3(d) are not hsted in the table because they are obtained from 
3(b) by permutations. They are the invariant eigenvectors of the same A4 members 
expressed in a similarly permuted basis. 

4. The first four solutions (3(b)) have the same energy, and the corresponding operators 
are all in class C3 or C4 of Table 2. The next three solutions (3(e)) also have the same 
energy, and the corresponding operators are all in class €2- In other words, vacuum 
alignments for operators in the same class always have the same energy. 

VI. GROUP THEORY AND DYNAMICS 

We saw in the last section that, barring permutations that serve only to relabel the 
neutrino mass eigenstates, every minimizing solution of the scalar potential V is an invariant 
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eigenvector of some group element of A4. Moreover, conjugate group elements give rise to 
identical energy. 

We shall present another derivation of these facts in the first subsection below. This 
alternative derivation relies more on the general property of groups, and is therefore more 
amenable to generalization. The generalization will be carried out in the second subsection. 

The following simple fact will be used in this alternative approach. 

Lemma: Let g"' be the a irreducible representation oi g E Let 0" be an eigenvector 
of g"" with eigenvalue A". Let :— (^"^^"^ • • • 0°")'' denote a state in the irreducible 

representation p obtained by coupling the n states 0"^. Then either = 0, or it is an 
eigenvector of g'' with eigenvalue A = 11^=1 -^"^ • 

Proof: This follows immediately from the fact that g'^{U.i (p'^'Y = {U.i{9°''(f>°'')y- 
The following is an immediate corollary of Lemma 1. 

A. Another derivation of the solution of U(l) x A4 

An alternative derivation of the result of the last section is presented here. It has the 
advantage of being easily generalizable from A4 to a general group Q. 

The main idea is the following. Start as before from the Higgs potential V — 
J2a fa{i^i^\'P°'\ipip) — fi"^ . Its minimization is determined by six vectorial equations 
of motion, dV/dipc = and dV/dip* = 0. The last three is the complex conjugate of the 
first three so they need not be considered separately. We will show, using the Lemma proven 
above, that if ijj is an invariant eigenvector of an element of A4, then every non-zero term 
in the three equations of motion is proportional to ■0, with a scalar coefficient common to 
the three equations. ip is then a solution of the equations of motion if we set this scalar 
coefficient to be zero, a condition that would be used to determined the normalization of 
the solution tjj. Once ijj is known, its energy can be computed from the above expression of 
V. 

Here are the details. The idea is simple but unfortunately the notations are a bit cum- 
bersome. 

has four irreducible representations, 1, 1', 1", 3. Other than the antisymmetric cou- 
pling 3 X 3 — )■ 3^ which does not appear in (31), and will therefore be forgotten from now 
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on, each irreducible representation a occurs at most once in the decomposition of ^ x 7. 

The Clebsch-Gordan coefficients (a a|/36, 7c) for /3 x 7 — >■ a can be expressed in terms of 

(a /3 7\ 

the 3j-symbol, {a a|/36, 7c) = \ [a] , where [a] is the dimension of the irreducible 

^ \a b c) 

representation a, and the 3j-symbol is completely symmetric in its columns. The state \oi a) 
is the complex-conjugated state of \aa). See Appendix A for further details. 

Consider a typical quartic term in the real Higgs potential V — J2a fa{'^'^l^\'P°'\i'fp) — 

(3 3 Oi\* { ol 3 3\ 
\\i>^)\i>o). (40) 
e/a/ \a b cj 

Its contribution to the equation of motion is 

Bn°' /3 3 a\* /a 3 3\ 

and its complex conjugate dQ" /dtp*. In these formulas, 1-0;,) stands for the 6th component 
of a three-dimensional complex column vector ip, and ((ipe\ the complex conjugate of \il)e)- 

We will use gf to denote IRp of a e A4. Suppose g"" has an orthonormal set of 
eigenvectors 0°^"^ (m = 1, 2, • • • , [a]), with eigenvalues A^, then 5 ad = Em 4>T^*4'd^^ so (40) 
and (41) can be rewritten as 

Q° = [a]^r"'"*F"'", (42) 

TO 

<5c = 2[a] ^ 2[q;] ^ (43) 

m TO 

where 

/a 3 3\ 

y"- = ^ Oo^c, (44) 

abc \a b c) 
(a 3 3\ 

= E (45) 

aft \a 6 c/ 

It follows from the symmetry of the 3j-symbol that y""* transforms as 1, and transforms 
as 3. 

We can now start working towards the goal of showing that invariant eigenvectors of g 
are solutions to the equations of motion derived from V . 

Let ip be an invariant eigenvector of g^ i.e., g^t/j — ip. Since gr"0""^ = A^0""^, according 
to the Lemma, either y""^ = 0, or gr^y""^ — A^F""*. However, since g^ — 1 for every g in 
the identity representation 1, we must have Aj„ = 1 in order for y""* 7^ 0. 
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The Lemma also implies that either ^""^ = 0, or gr^^am — A^^""*- Consequently, if 
Qam _^ meaning both y""* and non-zero, then it is necessary for to be an 
invariant eigenvector of g^. We know from Table 8 that every has one and only one 
invariant eigenvector, hence = K'^ip for some 7^ 0, if Q"™ 7^ 0. 

The importance of this observation is the following. There are six equations of motion, 
namely, the following three [c— 1, 2, 3), and their three complex conjugates: 

«=11'1"3 rn a=ll'l"3 "t 

From the observation above, we see that if ■0 is an invariant eigenvector of g^, then every 
non-zero term in (46) is proportional to ■0, hence all the six equations of motion are satisfied 
provided the single scalar equation 

E 2[a]/,F"™*«:"-/.2 = (47) 

a=ll'l"3 

is obeyed. Note that the single m in the equation is that one for which is the invariant 
eigenvector of g^. Since every quantity in this equation is known, except the overall normal- 
ization of ■?/', this then is an equation to determine the normalization of ip. For that reason 
we shall sometimes refer to (47) as the normalizafion equation. 

Knowing that -0 is an invariant eigenvector of g with known normalization computed in 
(47) , we can substitute this solution into the expression of V to compute its energy. 

To summarize, we have shown that the invariant eigenvector of every group element g 
of Ai is a dynamical vacuum alignment, satisfying the equation of motion derived from the 
Higgs potential V . We will now show that conjugate elements in give rise to the same 
energy V . 

Let g' = hgh~^ be an element conjugate to g. If g^ijj — -0, then g'^tp' — "0' for i/j' — h^ijj. 
Similarly, ^'«(/,'""» = A^^'""' if ^"0'*"* = A^^""' and 0'^*"* = /i°0""^. If we denote the 
quantities in (44) and (45) built up of the primed quantities by y'""* and ^'^^ respectively, 
then Y'""^ = and = /i^^""'. Furthermore, if = k"^, then ^'""^ = K^'ip' , with 
the same k", For that reason, the normalization equation (47) is identical for the primed as 
the unprimed quantities, so ifj and if]' have the same normalization, thus both solutions give 
rise to the same energy V . 

The eigenvalues of are — 1, — 1, and the eigenvalues of all other gf can be obtained 
from Table 8. Using this information, we can make a closer comparison between the present 
derivation and the explicit calculation in Sec. V. 
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If g is an order-3 element, then g^" 7^ 1 in the a — 1', 1" representations, hence Y°''^ — 
for a — 1', 1". This corresponds to the yy — yv — solutions in the previous section. 
Namely, the first four rows of Table 10, or solution 3(b) before that table. 

If g is order-2, for the singlet representations a = 1, 1', 1", every = +1, hence according 
to the Lemma, none of these three y""* (there is only one m) have to vanish, and we can 
take ^""^ = 1 in all cases. This corresponds to the solutions with \yi\ — \yy \ — \yv'\ 7^ in 
the last section, given by 3(e) and the next three rows of Table 10. 

B. Generalization to other groups Q 

The alternative derivation for discussed above is so general that it can be easily 
adapted to many other groups Q. For technical simplicity we shall confine ourselves here to 
the simply reducible groups. The validity of the theorem below certainly goes beyond the 
simply reducible groups, as A4 itself is not simply reducible. For those groups, the same 
idea should still be applicable, though generalization and some modification in details may 
be necessary. 

A simply reducible group is a group in which (i) g and g~^ are always in the same class, 
and (ii) in the decomposition of the Kronecker product /3 x 7 of any two irreducible repre- 
sentations /3 and 7, every irreducible representation a occurs at most once. A consequence 
of (i) is that all characters have to be real. A4 is not simply reducible because it contains 
complex characters, but groups like SU{2), SO{3), S3, S4, and the quaternion group, are 
simply reducible. 

The main property we need from a simply reducible group is the existence of a symmetric 
or antisymmetric 3j-symbol. As before, if {aa\f3b, 7c) is the Clebsch-Gordan coefficient, then 
the 3j-symbol is defined by 



is either symmetric or antisymmetric upon the interchange of columns. Here (a a) is the 
complex- conjugate representation of (aa); whether the 3j-symbol is symmetric or antisym- 
metric is a detail that does not concern us. 

Most of the arguments used in can now be carried through to prove 
Theorem: Let be a complex variable vector in the p irreducible representation of a 




(48) 
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simply reducible group Q, carrying a non-trivial U{1) quantum number opposite to that of 
Writing the most general {U{1) x ^)-invariant renormahzable real Higgs potential as 

V = ^/^(V'V'I^IV'V') - l^'im, (49) 

a 

in which the sum is carried over all irreducible representations, with being the projection 
operator and fa the Yukawa coupling constant of the a representation. Then every unique 
invariant eigenvector ip of g'' for every g & Q isa solution of the equation of motion dV/ dipc — 
dV/dtjj* — 0. Moreover, gr's in the same conjugacy class have the same Higgs energy V. 

Proof: The proof is essentially a copy of the procedure discussed in the last subsection. 
The only thing we are not sure without knowing more details about Q is whether every g^ 
has one and only one eigenvector with eigenvalue +1. This is why the theorem is limited 
to those gP with unique invariant eigenvectors. The only other difference is that the 3j- 
symbol may be antisymmetric under a column exchange rather than symmetric. The only 
place this symmetry is used was in (44) and (45) to show that they transform as 1 and 
3 respectively. Clearly an additional minus sign under column exchange will not affect 
this property. Finally, whether i}) belongs to 3 or some other multi-dimensional irreducible 
representation p really does not matter in the proof. 

With these preliminary remarks, let us proceed with the details. 

Consider a typical quartic term in the real Higgs potential V , 

^ {^^\V^\^^)^[a] (V'eKV'/lf'' r ^ ^\^^)\^c). (50) 

a,b,c,ej \e / a J \a b cj 

Its contribution to the equation of motion is given by 

80°' I P P p p\ 



^■^c ^ a,fe,c,e,f \e f a) \a b 



c , 



together with its complex conjugate dQ"/dip*. 

Suppose g"' has an orthonormal set of eigenvectors 0""* (m = 1, 2, • • • , [a]), with eigen- 
values A^. Then Sad — J2m 4>a^* 4>d^ : SO (50) and (51) can be rewritten as 

Q" = [q;]^^""**^"'", (52) 
m 

= 2[a]$:y"-*$r 2[a]Egr, (53) 
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where 



Y' 



am 




(55) 



(54) 



It follows from the symmetry of the 3j-symbol that y"*" transforms as 1, and transforms 
as p. 

Suppose tl) is an invariant eigenvector of i.e., g^ip — ip, and suppose this is the only 
eigenvector of gf with eigenvalue +1. Since g°'(f)°''^ — A^^""*, according to the Lemma, 
either = 0, or g^Y°''^ = XmY'^"^. However, since g^ ~ 1 for every g in the identity 
representation 1, we must have Am = 1 in order for ^ 0. 

The Lemma also implies that cither = 0, or gp^"'"^ = Am$"™. Consequently, if 

Qam ^ 0, it is necessary for to be an invariant eigenvector of g^^. Since this is assumed 
to be unique, it follows that = ^"^■0 foi" some k° if Q""* 0- 

Let us look at the equation of motion: 



From the observation above, we see that if ip is an invariant eigenvector of g/^, then every 
non-zero term in (56) is proportional to ip, hence equation of motion are satisfied provided 
the single scalar equation 



is obeyed. Note that the single m in the equation is that one for which is the invariant 
eigenvector of g''. Since every quantity in this equation is known, except the overall normal- 
ization of ijj, this then is an equation to determine the normalization of ■0. For that reason 
we shall sometimes refer to (57) as the normalization equation. 

This proves the first part of the theorem. What remains is to show that conjugate 
elements in Q give rise to the same energy V. 

Let g' = hgh~^ be an element conjugate to g. If gf^i/j — -0, then g'^ip' — 0' for 0' — h'^ijj. 



Similarly, ^'«0'°"^ = A„0'""* if ^"0°^"^ = A^0°"* and 0'°^"^ = /i°^0°"*. If we denote the 
quantities in (54) and (55) built up of the primed quantities by y'""* and respectively, 
then y'""* = Y""^ and ^'""^ = hP^"'"'. Furthermore, if ^""^ = K°'ip, then = K°'ip', with 



(56) 



am am 



E2N/ay""**««-/x^ = o 



(57) 



a 
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the same k", For that reason, the normahzation equation (47) is identical for the primed as 
the unprimed quantities, so ^0 and ■0' have the same normahzation, thus both solutions give 
rise to the same energy V. 

VII. CONCLUSION 

We have discussed in some detail the group-theoretical connection between a neutrino 
mixing matrix and the horizontal symmetry of left-handed leptons. Horizontal symmetry 
groups can be derived from mixing matrices, and mixing matrices can be obtained from 
horizontal groups, all without invoking the Higgs or other dynamical mechanisms, as long as 
we assume some trace of the original symmetry remains to be found in the mass matrices even 
after symmetry breaking. We have also shown that this formulation is completely consistent 
with the dynamical approach, in that the mixing matrices derived from a symmetry group Q 
can be obtained from a generic Higgs potential invariant under C/(l) x ^, for many groups Q. 
Such a potential gives rise to vacuum alignments that are invariant eigenvectors of the group 
elements. An alignment should be assigned to the charged-lepton sector if the corresponding 
group clement has order > 3, but to the neutrino sector if the corresponding group element 
has order 2. Within this rule, alignments can be further discriminated by their energies. 

Other than Sec. IV, in which some 5*4 and dynamical models arc constructed to 
illustrate the consistency of the two approaches, we have not studied in detail in this article 
how physical models should be constructed. The models constructed in Sec. IV may not be 
optimal, in that mass hierarchy is not natural in those models, and the successful Koide mass 
relation of the charged-leptons [9] is not explained. More complicated dynamics including 
the presence of quadratic valon fields may be necessary. Furthermore, although the U{1) xQ 
potential can give rise to all the desired vacuum alignments, and although their assignments 
to the two sectors can be discriminated by the order of the associated residual symmetry 
operator, we have not found a purely dynamical way to do so without invoking the residual 
symmetry operators explicitly. That may call for the introduction of driver valons [10] or 
other dynamical mechanisms. 
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Appendix A 



In this appendix, we list the Clebsch- Gordon (CG) series of A4 and its CG coefficients 
computed in the representation of Table 8. 

The CG series ^ x 7 — > a can be computed from Table 2 to be 
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I" 
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1" 


1" 
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1' 
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1, 1', 1", 3a 



Table 11. Clebsch-Gordan series of 

where the first row refers to IR7, the first column refers to IR/3, and the table entries refer 
to IRa. The CG coefficients (a a|/36, 7c) are zero if (a, /3,7) is not contained in an entry 
of the table. It is also zero if the product of F-eigenvalues of \^h) and |7c) is not equal to 
that of I a a). This last condition is a consequence of eq. (22) applied to g — F. If we use 
A, B, C to denote the three states in 1, 1', 1", and D, E, G to denote the three states in 3, 
then the F-eigenvaluc of A and D is 1, that of B and E is cj, and that of C and G is u"^. The 
eigenvalue rule is not shown in Table 11, but that is incorporated in the expanded Table 12. 
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B, Es 



Table 12. Expanded Clebsch-Gordan series of A4 

The complex conjugate \a a) of a state \aa) is defined to be the state with the complex- 
conjugated F-eigenvalue. Thus X = X if X = A, D, but S = C, C ^ B,E ^G,G ^ E. 

We will also write the CG coefficient simply as {X\Y, Z) to represent the more complicated 
{a a\f3b, 7c). 

41 



{a a 1/36, 7c) = a/ [a] 



The CG coefficients are related to tlie 3j-symbols via tfie equation 

a b c J 

{X\Y,Z) = Y Z), (Al) 

wfiere [a] is tfie dimension of IRa, and [X] is tfie dimension of tfie IR in whicfi state X 
belongs to. The 3j-symbol couples three states to a singlet; it vanishes unless the product 
of the three F-eigenvalues is 1. 

There are two 3j -symbols when a — f3 — ^ — 3. The one corresponding to a = 3^ will 
be marked with a subscript A, and the one corresponding to a = 35 will carry no subscript. 

3j-symbols are more convenient because of their symmetry properties. The 3j-symbol 
{X Y Z)a is antisymmetric upon the interchange of two columns, and all the other 3j- 
symbols are symmetric upon such an exchange. 

The non-zero 3j-symbol (X Y Z) is equal to 1 when X, Y. Z all belong to 1-dimensional 
IRs, is equal to l/y/3 if two of them belong to the 3-dimensional IR and the third one 
belongs to a 1-dimensional IR. Finally, when all three of them belong to the 3-dimensional 
IR, then their values are 

{D D D) = {E E E)^{GGG)^ (A2) 
It can be checked that the CG coefficients obey eq. (22) for g — G2 bs well. 
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